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Abstract. In this paper we compute the limit distributions of the numbers of hairpin-loops, 
interior-loops and bulges in fc-noncrossing RNA structures. The latter are coarse grained RNA 
structures allowing for cross-serial interactions, subject to the constraint that there are at most 
k — 1 mutually crossing arcs in the diagram representation of the molecule. We prove central 
limit theorems by means of studying the corresponding bivariate generating functions. These 
generating functions are obtained by symbolic inflation of lv|-shapes 



1. Introduction 



An RNA molecule is a sequence of the four nucleotides A, G, U, C together with the Watson-Crick 
(A-U, G-C) and U-G base pairing rules. The sequence of bases is called the primary structure 
of the RNA molecule. Two bases in the primary structure which are not adjacent may form 
hydrogen bonds following the Watson-Crick base pairing rules. Three decades ago Waterman et 
al. [HIES EE] analyzed RNA secondary structures. Secondary structures are coarse grained RNA 
contact structures. They can be represented as diagrams and planar graphs, see Fig. [TJ Diagrams 
are labeled graphs over the vertex set [n] = {1, . . . , n} with vertex degrees < 1, represented by 
drawing its vertices on a horizontal line and its arcs (i < j), in the upper half-plane, see 

Fig. [T]and Fig. [21 Here, vertices and arcs correspond to the nucleotides A, G, U, C and Watson- 
Crick (A-U, G-C) and (U-G) base pairs, respectively. In a diagram two arcs and (^2,^2) 
are called crossing if ii < i 2 < ji < 32 holds. Accordingly, a fc-crossing is a sequence of arcs 
(ii, ji), . . . , (ik,jk) such that ii < i 2 < ■ ■ ■ < ik < ji < h < ' • • < jk, see Fig. [2j We call diagrams 
containing at most (k — l)-crossings, fc-noncrossing diagrams (fc-noncrossing partial matchings). 
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FIGURE 1. The Sprinzl tRNA RD7550 secondary structure represented as 2-noncrossing 
diagram (top) and planar graph (bottom). 

The length of an arc is given by j — i, characterizing the minimal length of a hairpin loop. A 
stack of length r is a sequence of "parallel" arcs of the form 

(1-1) (i + 1, j - 1), (i + (r- 1), j - (r - 1))), 

and we denote it by SJj. We call an arc of length one a 1-arc. A fc-noncrossing, r-canonical RNA 
structure is a fc-noncrossing diagram without 1-arcs, having a minimum stack-size of r, see Fig. [5J 
Let 7k,T {n) denote the set of fc-noncrossing, r-canonical RNA structures of length n and let Tk,r {n) 
denote their number. 

We next introduce the following structural elements of fc-noncrossing, r-canonical RNA structures, 
see Fig. [3] and Fig. 2J Let denote an interval, i.e. a sequence of consecutive isolated vertices 
(i,i + 1, . , . ,j — 1, j). We consider, see Fig. @] 

(1) a hairpin-loop is a pair 

(2) an interior-loop is a sequence 

[k + 1, *2 - 1], (»2,j2), \h + l,ii - 1]), 
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FIGURE 2. A 2-noncrossing, 2-canonical RNA structure (left) and a 3-noncrossing, 2- 
canonical RNA structure (right) represented as planer graphs (top) and diagrams (bot- 
tom). 





110 



FIGURE 3. 3-noncrossing, 6-canonical structures: the pseudoknot structure of the PrP- 
encoding mRNA represented as diagrams (top) and planer graphs (bottom).. 



where («2, J2) is nested in (ii, ji). 

(3) a bulge is a sequence 

[ii + M2 - 1], («2,ji - 1)) or ((ii, ji), (ii + l,j 2 ), b*2 + l = - !])■ 

(4) a siem is a sequence of stacks 
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where the stack SJ m „■ is nested in Sj m 1 „■ , 2 < m < s and there are no arcs of the 

tjn ,Jm 'm — 1 , Jm — 1 

form (ii - 1, ji + 1) and (i s + r s , j s - t s ). 




1 23456789 10 11 12 13 




FIGURE 4. The loop-types: hairpin-loop (top), interior-loop (middle) and bulge (bottom). 

2. Preliminaries 

Let fk{n,£) denote the number of fc-noncrossing diagrams on n vertices having exactly I isolated 
vertices. A diagram without isolated points is called a matching. The exponential generating 
function of fc-noncrossing matchings satisfies the following identity [2J IH E] 

2n 

(2-1) H fe (z) = A(2»,0) • T^ry = det[/ w (2z) - /i +J -(2«)]|*ji 1 

n>0 ^ 

where I r (2z) = J2j>a j\(j+r)\ * s * ne hyperbolic Bessel function of the first kind of order r. Eq. (|2.ip 
allows us to conclude that the ordinary generating function 

n>0 

is D- finite [T^]. This follows from the fact that I r (2z) is D-finite and D-finite power series form 
an algebra [T2]. Consequently, there exists some e £ N such that 

(2.2) q Q>k {z)^F k {z) + g 1 , fe (z)^ rT F fc (z) + • • • + q e>k {z)¥ k {z) = 0, 

where qj >k (z) arc polynomials and qo. k (z) ^ 0. The ordinary differential equations (ODE) for F k (z), 
where 2 < fc < 7 are obtained by the MAPLE package GFUN from the exact data of f k (2n, 0). They 
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k qo,k(z) Rk 

2 (42-1)2 {1} 

3 (16z-l)z 2 {±} 

4 (144^-402 + 1)^ {±,^} 

5 (10242 2 - 80z + l)z 4 {je^} 



6 (14400z 3 - 4144z 2 + 1402 - l)z 5 

7 (1474562 3 - 125442 2 + 2242 - l)z 6 {ii^} 



TABLE 1. We present the polynomials qo,k(z) and their nonzero roots obtained by the 
MAPLE package GFUN. 

are verified by first deriving the corresponding P-recursions [12] for f^(2n, 0) second transforming 
these P-recursions into P-recursions of A-(2n, 0)/(2n)! and third deriving the corresponding ODEs 
for Hfe(z) and verifying that the RHS of eq. (|2.1[) is a solution. The key point is that any singularity 
of Ffe(z) is contained in the set of roots of qo,k(z) [H], which we denote by Rk- For 2 < k < 7, we 
give the polynomials qo,k{z) and their roots in Table [1] In [8] we showed that for arbitrary k 

(2.3) / fc (2n, 0) ~ c k n -((^-D 2 +(fc-D/2) ( 2 (fc _ !))*» ^ > 

in accordance with the fact that Ffc(^) has the unique dominant singularity p\, where p k = I /(2k — 
2). 

We next introduce a central limit theorem due to Bender pQ. It is proved by analyzing the 
characteristic function by the Levy-Cramer Theorem (Theorem IX. 4 in [3]). 

Theorem 1. Suppose we are given the bivariate generating function 

(2.4) f(z,u)= J2 f{n,t)z n u t , 

n,t>0 

where f(n,t) > and f(n) — ^2 t f(n,t). Let X„ be a r.v. such that P(X n = t) = f(n,t)/f(n). 
Suppose 

(2.5) [z"]/(z,e s )^ C ( S )^ 7 ( S )-" 

uniformly in s in a neighborhood o/O, where c(s) is continuous and nonzero nearO, a is a constant, 
andj(s) is analytic near 0. Then there exists a pair (/x, a) such that the normalized random variable 



6 



MARKUS E. NEBEL f , CHRISTIAN M. REIDYS* AND RITA R. WANG * 



has asymptotically normal distribution with parameter (0, 1). That is we have 

1 f x 

(2.7) lim P(X; < x) = —= \ e~i c2 dc 

™^°° V27T J-oo 

where fi and a 1 are given by 

(2.8) md V ' ( °> 



7(0) V 7(0) ) 7(0) ' 

The crucial points for applying Theorem [T] are (a) eq. (|2.5p 

[z n ]f(z,e s )~c(s)n a 7 (s)- n , 

uniformly in s in a neighborhood of 0, where c(s) is continuous and nonzero near and a is a 
constant and (b) the analyticity of 7(5) in s near 0. In the following, we have generating functions 
of the form Fk{ip{z, s)). In this situation, Theorem [5] below guarantees under specific conditions 



continuous, 



[z n ]FM(z,s)) ~ A(s) n-« fe - 1)2+ ( fe - 1 ^ 2 ' (-^] , A(s) 

for 2 < k < 7. The analyticity of j(s) is guaranteed by the analytic implicit function theorem [3]. 
Theorem 2. [7] Suppose 2 < k < 7 . Let ip(z, s) be an analytic function in a domain 

(2.9) T> = {(z,s)\\z\ <r,\s\ < e} 

suc/i f/ia£ "0(0, s) = 0. /n addition suppose j(s) is the unique dominant singularity of Fk(ip(z, s)) 
and analytic solution of i^(j(s), s) = pi, \j(s)\ < r, d z tj}{'~f{s), s) 7^ for \s\ < e. Then F%{tp(z, s)) 
has a singular expansion and 

(2.10) [z n ]FM(z, s)) ~ A(s) n -((fe~i) 2 +(fc-i)/2) ( J_\ f or some con Unuous A(s) £ C, 

\7(s)y 

uniformly in s contained in a small neighborhood ofO. 



To keep the paper selfcontained we give a direct proof of Theorem [5] in Section [5j This avoids 
calling upon generic results, such as the uniformity Lemma of singularity analysis [3]. 
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3. The generating function 

In this section we compute the bivariate generating functions of hairpin-loops, interior- loops and 
bulges. Let hk, T {n,t), i& )T (n,i) and bk >T (n,t) denote the numbers of fc-noncrossing, T-canonical 
RNA structures of length n with t hairpin-loops, interior-loops and bulges. We set 

(3.1) H fc , T («,ui) = ^^/i fe , T (n,t>"< 

n>0 t>0 

(3.2) Ife,r(2,u 2 ) = ^^i fe , r (n,t)2 n 4, 

ra>0 t>0 

(3.3) B fe , T (z,u 3 ) = ^^& fc ,.(MK l u 3 . 

ra>0 i>0 

In order to derive the above generating functions we use symbolic enumeration [3] . A combinatorial 
class is a set of finite size with the definition of size function of its elements, whose elements are 
all finite size and the number of certain size elements is finite. Suppose C be a combinatorial class 
and c e 6. We denote the size of c by \c\. There are two special combinatorial classes £ and Z 
which respectively contains only an element of size and an element of size 1. The subset of C 
which contains all the elements of size n in C is denoted by C„ . Then the generating function of a 
combinatorial class C is 

(3.4) C(z)=5>l c l = ]Tt7„z™, 

cee n>a 

where C n C C and C n = |C n |. In particular the generating functions of £ and Z are given by 
E(z) = 1 and Z(z) = z. For any two combinatorial classes 6, D, we have the following operations: 

• e + D:=euD, ifenD = 

• e x T> := {(c,d)\c ee,deV} and G rn := UT=i e 

• seq(c) = £ + e + e 2 + ---. 

We have the following relations between the operations of combinatorial classes and the operations 
of their generating functions: 

(3.5) A = e + T> => A(z) = C(z) +D(z) 

(3.6) A = £xT) => A(z) = C(z) • D(z) 

(3.7) yt = SEQ(C) A(z) = (l-C(z))- 1 , 
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where A(z), C(z), D(z) is the generating function of A, C and 3D. 

Given a fc-noncrossing, r-canonical RNA structure S, its lv|-shape, Ivf. (J) [TT], is obtained by first 
removing all isolated vertices and second collapsing any stack into a single arc, see Fig|3] By 
construction, Iv^-shapes do not preserve stack-lengths, interior loops and unpaired regions. In the 
following, we shall refer to Iv^-shape simply as shape. Let Tfc jT denote the set of fc-noncrossing, 
r-canonical structures and 3 k the set of all fc-noncrossing shapes and 3k(m) those having m 1-arcs, 
see Figure [5j Each stem of a fc-noncrossing, r-canonical RNA structure is mapped into an arc 
in its corresponding shape and all hairpin-loops are mapped into 1-arcs. Therefore we have the 
surjective map, 

(3.8) <p: 7 Kr ^3k. 

Indeed, for a given shape 7 in 3k, we can derive a fc-noncrossing, r-canonical structure having 

1 2 3 4 5 6 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 1617 18 19 20 



FIGURE 5. A 3-noncrossing, 2-canonical RNA structure (top-left) is mapped into its 
shape (top-right) in two steps. A stem (blue) is mapped into a single shape-arc (blue). 
A hairpin-loop (red) is mapped into a shape-l-arc (red). 



arc-lcngth> 2, we can add arcs to each arc contained in the shape such that every resulting stack 
has r arcs and insert one isolated vertex in each 1-arc. Let 3k(s,m) and ik(s,m) denote the set 
and number of the Iv^-shapcs of length 2s with m 1-arcs and 

s 

(3.9) i k {x lV ) = Y J Y. i ^ s ' m ) xS y m 

s>0 m=0 

be the bivariate generating function. Furthermore, let Ufc(m) denote the set of shapes 7 having m 
1-arcs. Let fc,s,m be natural numbers where fc > 2, then the generating function lk(x, y) |llj is 
given by 

1 + x ( x(l + x) 



^ ^ k ^ ' ^ 1 + 2x — xy K \ (1 + 2x — xy) 2 
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Theorem 3. Suppose k, r e N, k > 2, r > 1. TTien 

(l-z)(l-z 2 + z 2r ) 



(3.11) 



(3.12) 



(3.13) 



Hfc !r (z,Ml) 



(1 - z) 2 (l - z 2 + z 2 ^) + z 2 ^ - z^+iui 

F ( z 2 -(l-z) 2 (l-z 2 + z 2 -) 

V ((1 - z) 2 (l - z 2 + z 2t ) + z 2t - z^mf 



(1 - z 2 )(l - z) 2 - u 2 z 2t + 2 + (2z 2 - 2z + l)z 2r 
(1 - z) ((1 - z 2 )(l - z) 2 - u 2 z 2t + 2 + (2z 2 - 3z + 2)z 2t ) 

' z 2t ((1 - z 2 )(l - z) 2 - M2 z^+ 2 + (2z 2 - 2z + l)z 2 ") ' 
((1 - z 2 )(l - z) 2 - u 2 z 2t + 2 + (2z 2 - 3z + 2)z 2 ^) 2 



F, 



(l-z 2 )(l-z)-2^ 2r+1 + (^ + l)^ r 
fe ' rlZ ' " 3j (1 - z) ((1 - z 2 )(l - z) - 2u 3 2 2r+1 + (^ + 2)z 2 ^) 

/ z 2r ((l-z 2 )(l-z)-2u 3 z 2r+1 + (z + l)z 2T ) 

k \ (1 - z) ((1 - z 2 )(l - z) - 2u 3 z 2r+1 + (z + 2)z 2T ) 2 



Proof. We prove the theorem via symbolic enumeration representing a fc-noncrossing, r-canonical 
structure as the inflation of a shape, 7. Since a structure inflated from 7 £ 3k{s,m) has exactly s 
stems, (2s + 1) (possibly empty) intervals of isolated vertices and to nonempty such intervals we 
rewrite the generating functions as 

Hfc, T (z,Ui) = ]P ^2 T 7( Z > U 1>M)> 

m>0 76 Jfc(m) 

Ifc lT (^,W2) = zJ S T ^ Z ' 1 ''" 2 ' ^ 
m>0 76 Jfc(m) 

B fc , T (z,u 3 ) = X] X] Tt^' 1 ! 1 '^)- 

m>0 T £3 t (m) 

where T 7 (z, ui, U2, u 3 ) is the generating function of all fc-noncrossing, r-canonical structures with 
shape 7 and u, (i = 1, 2, 3) are variables associated with the number of hairpin-loops, interior-loops 
and bulges. In order to compute the latter we consider the inflation process: we inflate 7 £ 3fc(m) 
having s arcs, where s > to, to a structure as follows: 



• we inflate each arc of the shape to a stem of stacks of minimum size r. Any isolated vertices 
inserted during this first inflation step separate the added stacks, 

• we insert isolated vertices at the remaining (2s + 1) positions. 
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We inflate any shape-arc to a stack of size at least r and subsequently add additional stacks. The 
latter are called induced stacks and have to be separated by means of inserting isolated vertices, see 
Fig. [51 Note that during this first inflation step no intervals of isolated vertices, other than those 
necessary for separating the nested stacks are inserted. After the first inflation step we proceed 




1 2 3 4 5 6 7 8 9101112131415 1617 



1 2 3 4 5 6 1 2345678910 1112 123456789101112131415 1617 




(3) 



1 2 3 4 5 6 7 8 910111213141516 1718 



FIGURE 6. The first inflation step a shape (left) is inflated to a 3-noncrossing, 2- 
canonical structure. First, every arc in the shape is inflated to a stack of size at least two 
(middle), and then the shape is inflated to a new 3-noncrossing, 2-canonical structure 
(right) by adding one stack of size two. There are three ways to insert the isolated 
vertices. 

inflating further by inserting only additional isolated vertices at the remaining (2s + 1) positions 
in which such insertions are possible. For each 1-arc at least one such isolated vertex is necessarily 
inserted, see Fig. [JJ We proceed by expressing the above two inflations in terms of symbolic 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 



FIGURE 7. The second inflation step: the structure (left) obtained in (1) in Fig. [6] is 
inflated to a new 3-noncrossing, 2-canonical RNA structures (right) by adding isolated 
vertices (red). 

enumeration. For this purpose we introduce the combinatorial classes M (stems), 3C T (stacks), N T 
(induced stacks), £ (isolated vertices), 3^ (arcs) and Z (vertices), where Z(z) = z and R(z) = z 2 . 
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Let hi, fi2 and /13 be the labels for hairpin- loops, interior- loops and bulges, respectively. Then 



(3.14) 




= {M) s x Z 2s+1 - m 


(3.15) 


M 


= % T x Seq(?T), 


(3.16) 


W 


= % T x f[Z x £] M3 


(3.17) 


% T 


= 3T x Seq (3?) , 


(3.18) 


L 


= Seq (Z) . 



and consequently, translating the above relations into generating functions the generating function 
T 1 (z,ui,u 2 ,u 3 ) is given by 

\ 



1^ 



1 - z 



2s+l-m 



Ml X 
1 - 2 



(1 Z) 1 ( v (1-z 2 )(1-z) 2 -(2 U3 z(1-z) + U2^ 2 )z 2t J (" lZ ) m ' 
where the indeterminants itj (i = 1, 2, 3) correspond to the labels fii; i- e - the occurrences of hairpin- 
loops, interior-loops and bulges. Accordingly, for any two shapes 71,72 € Jfe(m) having s arcs, we 
have 



(3.19) 
We set 
(3.20) 



T 7i (z,u 1; m 2 ,u 3 ) = T 12 (z,u 1 ,u 2 ,u 3 ) 



(1 - z 2 )(l - z) 2 - (2u 3 z (1 - z) + u 2 z 2 )z 2 



and accordingly derive 



H fe , T (z,ui) = X] T 7 (z,ui, 1, 1) = ^ ^ i fc (s,m)T 7 (z,Mi, 1, 1), 

m>0 7 G Jk(m) 



s>0 m=0 



lk, T {z, u 2 ) = ^2 ^2 T 7 (z, l,u 2 , 1) = ^ ^2 ife(s,m)T 7 (z, l,u 2 , 1), 

m>0 7 e3fc(m) s>0m=0 

s 

B fe , T (z,u 3 ) = ^2 T t( z AAiU 3 ) = ^2^2i k (s,m)T 7 (z,l,l,u 3 ). 

m>0 7 SJ fc (m) s>0m=0 

It now remains to observe 

1 + x „ 



£ £ i fc (s,m)x s y™ = 



s>0 m=0 



x(l + x) 
1 + 2a; — xy V. (1 + 2x — xy) 
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and to subsequently substitute x = 77(1,1) and y = u\z for deriving Hfc iT (z,iti). Substituting 
x = r/(u2,l) and y = z in we obtain l)- jT (z, u-x) and finally x = 77(1,1*3) and y = z produce the 
expression for Bfc )T (z, U3), whence the theorem. □ 



4. The central limit theorem 



For fixed fc-noncrossing, r-canonical structure, S, let M n ^, T (S), ^n,k,r(S) and B„.fc. r (S) denote the 
number of hairpin-loops, interior-loops and bulges in S. Then we have the r.v.s 

• H n , fc , T , where P(H n>fejT = t) = 

• I„ jfejT , where P(I„, fc , T = f) = 

• l n , fc , T) where P (B„. fe . T = t) = 

Here hk. T (ji,t), ik.r{n,t) and bk. T (n,t) are the numbers of fc-noncrossing, r-canonical structures of 
length n with t hairpin-loops, interior-loops and bulges. The key for computing the distributions 
of the above r.v.s are the bivariate generating functions derived in Theorem [3] 

(4.1) H fc , T («,ui) = ^2^2hk, T (n,t)z n u\, 

n>0 t>0 

(4.2) I k)T (z,u 2 ) = ^2^2i k , T (n,t)z n u t 2 , 

n>Q t>0 

(4.3) B Kt (z,u 3 ) = ^^^(MK 1 ^. 

n>0 t>0 

The following proposition is based on Theorem [2] and facilitates the application of Theorem [TJ 

Proposition 1. Suppose 2 < k < 7 , 1 < t < 10. There exists a unique dominant Hfe T (z, e s )- 
singularity, 7fc jT (s), such that for \s\ < e, where e > 0: 

(1) 7fc, T (s) is analytic, 

(2) 7fc. T (s) is the solution of minimal modulus of 

, . z 2T (l-z) 2 (l -z 2 + z 2t ) 2 
4.4 i i-i i 3 - p? = 0. 



-r 2 

((1 - z) 2 (l -z 2 + z 2T ) + z 2T - z 2T + 1 e s 

and 



(4.5) [z"]H fc , T (z,e s )~C( S )7i- 



7fe,r(s) 
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uniformly in s in a neighborhood of and continuous C(s). 

Proof. The first step is to establish the existence and uniqueness of the dominant singularity 

7fc,r(s). 
We denote 

(4.6) #{z,s) = (l-z) 2 (l-z 2 + z 2T ) + z 2T -z 2r +V, 

(4.7) M*,s) = z 2t (1-z) 2 (1-z 2 + z 2t )#(z,s)- 2 , 

(4.8) u T (z,s) = (l-z)(f-z 2 + z 2T )i?(z,s)" 1 , 

and consider the equations 

(4.9) V2<i<fc; Fi, T (z,a) = il> T (z t s) - f%, 

where /O, = l/(2i — 2). Theorem [3] and Table [T] imply that the singularities of Hj,. T (z,e s ) are are 
contained in the set of roots of 

(4.10) F ljT (z,s) = and tf(z,s) = 
where i < k. Let r 2 ;. T denote the solution of minimal modulus of 

(4.11) #,t(2,0) = VvM)-p 2 =o. 

We next verify that, for sufficiently small > 0, |z — r< T | < |s| < ej, the following assertions 
hold 

• ^i,T(ri,„0) ^0 

• J^F i T (z, s) and -^F i T (z, s) arc continuous. 

The analytic implicit function theorem, guarantees the existence of a unique analytic function 
7i, T (s) such that, for \s\ < e ; , 

(4.12) F i)T (j i)T (s),s)=0 and 7 ijT (0) = r i>T . 

Analogously, we obtain the unique analytic function S(s) satisfying $(z,s) = and where (5(0) 
is the minimal solution of $(z, 0) = for \s\ < eg, for some eg > 0. We next verify that the 
unique dominant singularity of Hfe lT (z, 1) is the minimal positive solution r^^ of Fk, T (z,0) = 
and subsequently using an continuity argument. Therefore, for sufficiently small e where e < e, 
and e < ea, |s| < e, the module of 7i, T (s), i < k and <5(s) are all strictly larger than the modulus 
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[z"]H,, r (z,e s )~C( S )n- 



of 7/b jT (s). Consequently, 7fc, r (s) is the unique dominant singularity of T (z, e s ). 

Claim. There exists some continuous C(s) such that, uniformly in s, for s in a neighborhood of 

1 

,7fe,r(s) / 

To prove the Claim, let r be some positive real number such that r k ^ T < r < <5(0). For sufficiently 
small e > and \s\ < e, 

l7fc,r(s)| < r and \5(s)\ > r. 
Then ip T (z,s) and lj t {z,s) are all analytic in D = {(z, s)\\z\ < r, \s\ < e} and "0 T (O, s) = 0. Since 
7fe,r(s) is the unique dominant singularity of 

H k , T (z, e s ) = uj t (z, s) F k (ip T (z, s)), 

satisfying 

(4.13) V , r(7fe,T(s),s) = Pk and \lk,r(s)\<r, 

for |s| < e. For sufficiently small e > 0, J^F k r (z 7 s) is continuous and -g^F k _ T (r k _ T , 0) 7^ 0. Thus 
there exists some e > 0, such that for \s\ < e, ■§^F kT {^ kT {s)^s) 7^ 0. According to Theorem^ we 
therefore derive 

(4.14) [z n ]n ktT (z, e s ) ~ C(s) n-(( fe " 1 ) 2 +^ i ) f—L 

\7fe,r(s), 

uniformly in s in a neighborhood of and continuous C(s). □ 



After establishing the analogues of Proposition!]] for l k , T (z, u) and Bfe jT (z, u), see the Supplemental 
Materials, Theorem [T] implies the following central limit theorem for the distributions of hairpin- 
loops, interior-loops and bulges in fc-noncrossing structures. 

Theorem 4. Let k, r G N, 2 < k < 7, 1 < r < 10 and suppose the random variable X denotes 
either M n- k. T , In.fc.r orB„.fc. T . TTien i/iere exists a pair 

sitc/i i/iai i/ie normalized random variable X* /ias asymptotically normal distribution with parameter 
(0, 1), where /ife,r,X arl ^ °1c T x are ff* wew &2/ 

7^,x(0) 2 / 7^(0) \ 2 7^,x(0) 

( j W,T,X ~ ^(o)' CTw_ l^xTo)J : t!^(o)' 

where 7k, Tj x(s) represents the unique dominant singularity ofHk. T (z, e s ), I^. T (z, e s ), andBfc j7 -(z, e s ), 
respectively. 
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k = 


2 


k 


= 3 


k 


= 4 














Hk,r 




T 


= 1 


0.105573 


0.032260 


0.012013 


0.011202 


0.003715 


0.003641 


T 


= 2 


0.061281 


0.018116 


0.009845 


0.008879 


0.003734 


0.003602 


T 


= 3 


0.043900 


0.012752 


0.007966 


0.007060 


0.003200 


0.003060 


T 


= 4 


0.034477 


0.009896 


0.006680 


0.005854 


0.002757 


0.002622 






k = 


5 


k 


= 6 


k 


= 7 










Vk,r 


^k.T 






T 


= 1 


0.001626 


0.001612 


0.000855 


0.000852 


0.000505 


0.000504 


T 


= 2 


0.001897 


0.001864 


0.001123 


0.001111 


0.000731 


0.000726 


T 


= 3 


0.001693 


0.001655 


0.001035 


0.001021 


0.000692 


0.000686 


T 


= 4 


0.001486 


0.001448 


0.000922 


0.000907 


0.000624 


0.000618 



TABLE 2. Hairpin-loops: The central limit theorem for the numbers of hairpin-loops in 
fc-noncrossing, r-canonical structures. We list fx^ >T and &1 >T derived from eq. (14.151) . 







k 


= 2 


k 


= 3 


k 


= 4 






^k,r 








^k,r 




T 


= 1 


0.015403 


0.013916 


0.001185 


0.001176 


0.000264 


0.000264 


T 


= 2 


0.012959 


0.011395 


0.001823 


0.001793 


0.000603 


0.000599 


T 


= 3 


0.011075 


0.009570 


0.001878 


0.001837 


0.000693 


0.000688 


T 


= 4 


0.009682 


0.008261 


0.001803 


0.001755 


0.000700 


0.000693 






k 


= 5 


k 


= 6 


k 


= 7 










Vk,r 




Hk,r 




T 


= 1 


0.000090 


0.000090 


0.000039 


0.000039 


0.000019 


0.000019 


T 


= 2 


0.000275 


0.000274 


0.000149 


0.000149 


0.000090 


0.000090 


T 


= 3 


0.000343 


0.000341 


0.000198 


0.000198 


0.000126 


0.000126 


T 


= 4 


0.000359 


0.000357 


0.000214 


0.000213 


0.000140 


0.000140 



TABLE 3. Interior-loops: The central limit theorem for the numbers of interior-loops in 
fc-noncrossing,r-canonical structures. We list fik,r and a\ T derived from eq. (14.151 1. 



In Tables [U [3] and d] wc present the values of the pairs (Mfc.r.x, cr| tX ). 
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k 


= 2 


k 


= 3 


k 


= 4 














Pk,r 




T 


= 1 


0.049845 


0.042310 


0.008982 


0.008684 


0.003094 


0.003058 


T 


= 2 


0.025088 


0.021785 


0.005789 


0.005597 


0.002457 


0.002422 


T 


= 3 


0.015859 


0.013979 


0.003936 


0.003814 


0.001762 


0.001737 


T 


= 4 


0.011197 


0.009980 


0.002878 


0.002795 


0.001318 


0.001301 






k 


= 5 


k 


= 6 


k 


= 7 










Vk,r 


^k.T 


Pk,r 




T 


= 1 


0.001422 


0.001414 


0.000770 


0.000767 


0.000463 


0.000462 


T 


= 2 


0.001326 


0.001316 


0.000817 


0.000813 


0.000547 


0.000546 


T 


= 3 


0.000991 


0.000984 


0.000632 


0.000629 


0.000436 


0.000435 


T 


= 4 


0.000755 


0.000750 


0.000489 


0.000486 


0.000342 


0.000341 



TABLE 4. Bulges: The central limit theorems for the numbers of bulges in fc-noncrossing, 
r-canonical structures. We list p,^ >T and a\ T derived from eq. (|4.15[) . 



5. Proof of Theorem [2] 



Proof of Theorem [2l Wc consider the composite function Fk(tp(z, s)). In view of [z n ]f(z, s) = 
7™[z n ]/(^-, s) it suffices to analyze the function Fk{ip{~f(s)z, s)) and to subsequently rescale in order 
to obtain the correct exponential factor. For this purpose we set 

where ip(z, s) is analytic in a domain D = {(z, s)\\z\ < r, \s\ < e}. Consequently ip(z, s) is analytic 
in \z\ < r and |s| < e, for some 1 < r, < e < e, since it's a composition of two analytic functions 
in D. Taking its Taylor expansion at z = 1, 

(5.1) i>{z,s) = Y,Ms){l-zT, 

n>Q 

where ip n (s) is analytic in |s| < e. The singular expansion of Ffc(z), 2 < k < 7, for z — > pi, follows 
from the ODEs, see eq. (|2.2[) . and is given by 

p = f Pk{z - pi) + c' k {z ^)((*-D a +(fc-i)/a)-i log(z - pi) (1 + o(l)) 
1 Pk(z - Pi) + c' k (z - ^)((*-D a +(fe-i)/a)-i (i + (1)) 
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depending on whether k is odd or even and where Pk(z) are polynomials of degree < (k — l) 2 + 
(fc— l)/2— 1, c' k is some constant, and p k = l/2(fc— 1). By assumption, 7(s) is the unique analytic 
solution of , 0(7(s), s) = p 2 and by construction Fk(ip("f(s)z 7 s)) = Fk{ip(z 7 s)). In view of eq. (|5.1|) . 
we have for z — > 1 the expansion 

(5.3) ${z, s) -pl = J2 - zf = Ms)(l - z)(l + o(l)), 

n>l 

that is uniform in s since ip n (s) is analytic for \s\ < e and ipo(s) = ^(l(s),s) = p\. As for the 
singular expansion of F k (ip(z, s)) we derive, substituting the eq. (|5.3[) into the singular expansion 
of F fc (z), for z-»l, 



(5.4) 



P fc (2, s) + c fc (s)(l - z )((fe-i) 2 +(fc-i)/2)-i i g(i - z) (1 + o(l)) for k odd 
P fe (z, s) + c fe (s)(l - z )((^i) 2 +(*-i)/2)-i (i + (i)) f or fc ev en 



where P k (z,s) = P k (ip(z,s) - p 2 k ) and c k (s) = c'^is)^-^ 2 ^- 1 ^- 1 and 

ipi(s) = d z ip(z,s)\ z =i = j(s)d z tp(j(s), s) ^ for |s| < e. 

Furthermore P k (z, s) is analytic at |z| < 1, whence [z™]Pfc(z, s) is exponentially small compared to 
1. Therefore we arrive at 



(5.5) [z n ]F k (r/,(z,8)) 



[z n }c k (s)(l - z )((fe-i) 2 +(fe-i)/2)-i i g(i - z )(l + o(l)) 
[z"]c fe (s)(l - z )(C=-i) 2 +(^-i)/2)-i (! + (!)) 



depending on k being odd or even and uniformly in \s\ <e. We observe that c k (s) is analytic in 
\s\ <e. Note that a dependency in the parameter s is only given in the coefficients c k (s), that are 
analytic in s. Standard transfer theorems [■'>■ imply that 

(5.6) [z n ]F k $(z, s)) ~ A(s) n -((fc-i) 2 +(fe-D/2) for somc ^( s ) G Cj 

uniformly in s contained in a small neighborhood of 0. Finally, as mention in the beginning of the 
proof, we use the scaling property of Taylor expansions in order to derive 

(5.7) [z n ]F k y,(z, s)) = ( 7 ( S ))- [z n )F k $(z, sj) 
and the proof of the Theorem is complete. 
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